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ABSTRACT 


Device Hiodeling has become an incispenssible tool 
in the design and development of complex LSI/VLSI circuits e 
The main object is to link the physical device parameters 
to the output terminal characteristics. An accurate model 
is used to predict rhe device behaviour without having to 
go into actual fabrication process. The overall result is 
the Considerable simplification in device design procedures 
and in the simulation of output performance. 

The packin:. density and performance level consiceratiq 
haveied to the design of k'.OSFBTs with very small dimensions. 
Such small geometry MOSFETs are used extensively in LSI/VLSI 
digital systems. These small geometry MOSrETscive rise to 
certain effects which complicate device operation and degradd 
device performance. It therefore becomes necessary and 
essential to understand the physical behaviour of such MOSFET 
and to develop suitable models. 



The modeling of higher order effects in small 
eometry MOSFETs, have been, thus carried out with above 
entioned aspects in view* The moceis developed to 
ccount the various higher order effects are based on 
emiempirical and analv'xical approach. The reason is to 
I'ovide a better insight into physical mechanism and also 
o develop faster computer program? 


suitable f or CAD 
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CHAPTER 1 


INTRCSXJCTION 


In the past several years MOS technology has become 
indispensable in the design and development of large scale 
integrated circuits and system. Using the MOS devices as 
building blocks, a vast array of complex integrated system 
has been designed. Because of its potential for high 
density, high performance and low power applications, MOS 
technology has been the tool for the rapid growth and 
advancement of MOS VLSI systems. 

In an effort to increase the silicon real estate 
utilization the integration density has beai increased 
steadily and when compared to the medium scale integration 
of a decade ago dramatically, significant part of this 
steady increase lies in the reduction of the size of the 
individual devices, and this reduction has been supported 
by several technological developments such as more accurate 
process control and fine pattern lithography by optical, 
electron beam, and ion beam techniques. However as the 
channel length is reduced, many effects, which, heretofore 
were of second order importance, became of primary importanc< 
and dominate device and circuit perfornaance. The reduction 
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of the device size in order to achieve greater performance 
has followed a scaling principle [ij* hut this approach is 
limited by physical and practical problems. An insight 
into some of these problems can often be obtained through 
the study of models of the device itself. 

The technological difficulties involved and the 

possible innovative solutions have been discussed by a 

number of authors [2] - [4], For the present case, main 

emphasis here is on the study and analysis of various small 

geometry MOS effects necessary to predict accurately 

overall device behaviour. One of the basic limitations of 

size reduction in VLSI is the spacing of the source and 

drain diffusions i,e, the channel length. Drain potential 

reverse biases the drain- substrate p-n junction and creates 

a field pattern that can lower the brarrier, separating it 

is 

from the source. As the barrier flowered, additional 
subthreshold current can flow, and the source-substrate-drain 
structure acts qualitatively as a bipolar transistor. For a 
given channel doping concentration, as the channel length 
is reduced, the depletion layer width of the source and the 
drain junctions become comparable to the channel length. The 
potential distribution in the channel now depends on both 
the transverse field E (controlled by gate and substrate 
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voltage) and longitudinal field (controlled by drain 
bias). This two dimensional potential distribution is 
in contrast with original gradual channel approximation 
(E >> E ) for long channel devices. This two dimensional 
potential results in degradation of various device para- 
meters. The present work emphasizes study of such higher 
order effects and to provide simple models based on 
reasonable assumptions. 

The study of higher order effects as such is important 
for the development of CAD models for the moSFET which would 
greatly help in more accurate prediction of the device and 
circuit performance. They would also aid in the simulation 
of scaled down MOSFETs. Besides it would also be possible 
to study the influence of specific device parameters on the 
output behaviour of the basic logic gates. Such studies 
Would significantly contribute towards optimum selection 
and design of the various device parameters. Finally since 
actual process dependent effects have been included in the 
models developed, some indirect link between the electrical 
behaviour and the fabrication process would also be 
obtained. 

In Chapter 2, a brief discussion on long channel 
MOSFET theory has been presented. The theory aims to provide 



a link between the existing long channel expressions with 
correspondiing short channel relationships described in later 
chapters. A short introduction to general MOSFET fabrication 
process is given. This is required to appreciate the various 
process dependent effects. 

Chapter 3 covers the discussion on threshold voltage 
of small dimension MOSFETs. This includes the models of 
threshold voltage for short channel and narrow width devices. 
Also modification due to nonuniformly doped substrate 
structure has been included. 

In Chapter 4, the various small dimension 'effects 
have been described, which tends to limit the useful operating 
voltage range of small geometry MOSFETs. 

Chapter 5 covers the conduction behaviour of small 
geometry MOSFETs at below and above threshold voltages. 

In Chapter 6, the influence of threshold voltage 
variation due to scaling and influence of subthreshold current 
on circuit performance is discussed. 

Finally Chapter 7 presents the concluding remarks 
and discussions on several relevant aspects. Scope for 
future Work has also been indicated. 



CHAPTER 2 


MOSFET THEORY 


2.1 INTRODUCTIOJ^ 

The schematic diagram of a typical n channel MOSFET 
is shown in Fig, 2.1, the coordinate system has also been 
shown, the direction normal to GATE has been taken as x axis 
whereas y axis is ta^en along the channel in the direction 
of current flow, the same convention has been used throughout 
the present work. 

Two n-type regions are diffused into the p type 
substrate these regions form the drain and source contacts. 
The gate structure is essentially combined with the p type 
substrate to form a MOS structure. For a enhancement mode 
device, with no voltage applied to the gate, the source to 
drain electrodes correspond* to two p-n junctions connected 
back to back. The only current that can flow from source to 
drain in this case is, the leakage or subthreshold current. 

If the gate is biased positively, a aegative surface space 
charge appears at the semiconductor surface next to the 
oxide. For a sufficiently large forward bias, a n type 
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inversion layer forms at the surface* This inversion layer 
forms a narrow channel betv^een the source and drain contactSj 
and this inversion channel conducts current from the source 
to drain. The MOS structure modulates this current by 
varying the surface charge. The depth of the channel into 
the p-region is determined by the gate voltage .and the drain 
voltage, since it is the difference between voltage accross 
the region - V(y) , where V(y) is the surface potential 
at y, which determines the surface chaifge density at that 
point* The backgate bias or substrate bias also affects 
the channel conductance, 

2,2 THE SURFACE POTENTIAL 

Formation of the channel at the surface will occur 
at strong inversion that is when the minority carriers at 
the surface become equal to the majority carrier density at 
the bulk Ug = this occurs.' when 

With = a in (^). 

For grounded source and drain voltage > 0, surf a c 
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potential vary from source to drain, corresponding energy 
band diagram is shovm in Fig, 2.2, obviously 


^s 

*inv 


y=0 


20 f 


and <p 

inv 


y=L 


= 20 f + 


From this it is apparent that the surface potential 
increases from source to drain, and voltage across oxide 
layer and channel is a decreasing function of y, as y increases 
from source to drain. The channel depth decreases along 
same direction with pinch off occurres at the drain end of 
the channel, 

2.3 THE THRESHOLD VOLTAGE 

The voltage applied across the gate appears partly 
across the oxide and partly across the semiconductor. Thus 


Vgs = 


0 

^MS ^s ” 


MS + ^ox ^^s 


ox 


= ( 0 ; 


MS 


*^ox 




( 2 . 1 ) 




Fig. 2.3 Formation ^Depletion Layer and Cha.nnl 
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where. 



= metal semiconductor work function difference 


Qp = Fixed oxide charges density 

Q- = Interface state charge density 
mL 3 

Qg = depletion layer charge density 
C „ = Oxide capacitance per unit area 

OX 


Threshold voltage, defined as gate voltage at onset of 
inversion with = 2 0^ 


V 


th 



+ 20 



( 2 . 2 ) 


where Vpg = flat band voltage, 

2.4 CURRENT VOLTAGE CHARACTERISTICS 

Assuming an NMOS device |b voltages V^g, V^g and Vgg 
applied as shown in Fig. S.S^jjjue to presence of Vds- ® 
wider depletion layer exists at drain. 

At a distance y,V(y) and gate to channel voltage at 
that point is Vq - V(y). If this voltage exceeds threshold 
voltage conducting channel is formed. The induced charge per 
unit area in the channel is at the point y. 


Qi(y) = 9sc^y) - o^Cy) 


(2.3) 
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where, 

= space charge layer charge per unit area, 

«SC = -<''g - Cox 

Thus Q.(y) = -(Vgs - %) C^y, + r(2q 

Now surface band bending at a distance y is given as 

% = % + V(y) 

®inv 

Thus Q, = -[Vq 3 - I V(y) + 3 

+ r 2q < (V(y)N +9s ) 

^ * * ®inv 

(2.4) 

Assuming the conducting channel, to be homogeneous, 
so that ohms law is valid. The distance dR of a length 
dy of the channel 

= W.lin|ai^)| 

where ^n = average mobrility of the electrons in the channel 


W = channel width 
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The voltage drop along the length of the channel dy is 
given by 

dV = Ij^dR 

Substituting Qj^(y) from (2.4) and integrating left side in 
voltage from V = 0 to V = and Right side in length from 

y = 0 to y = L, 


''DS In L , , 

f = riS- [ i =ox 


[2q6^N. 


inv 


'^n oxJu r 
= -T — rDsl 


" L 


^ BSL Vgg I ^ t(VQ3+ 


with k = 


V(2iMA' 


- ( Si’s ] ] 

®inv J 

(2.6) 


ox 


For small V^g, Eq, (2.6) is simplified as according 
to Shichman hedges model [5]: 



14 


“ L ^ox 


DS " 2 i 


(2.7) 


^DS 1 ^GS - ^th ^ 0 


where V+,. = S'*. + 

®inv 




FB 


OX 


for linear region of operation. 


And, 


Id = It^n^ox tVQS-Vth]^ 


( 2 . 8 ) 


"DS 


^GS - ^th i 0 


2.5 FABRICATION PROCESS 


The state of the art process sequence for fabrication 
of n channel AK)S^ET may be summarized in following steps [6], 

a) A chemical vapour deposition (CVD) process deposits a 
thin layer of silicon nitride (Si2N^) on the entire wafer 
surface, silicon nitride is removed selectively and Boron is 
implanted . in exposed areas to suppress unwanted conduction. 
Next a layer approx. Ip thick of silicon dioxide (Si02) is deposa 
in these inactive or field regions. 

b) The Si^N^ is removed and a clean thermal oxide about 0.1 
micron thick is grown. Another CVD process deposits a layer 
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of polysilicon over the entire wafer. Next photolichocra phy 
ctep defineS' the desired patterns for gate electrodes. Un- 
cesired poly is removed by chemical or plasma etching. An 
n x.ype dopant is introduced to form source or drain by 
either thermal diffusion or ion-implantation . process. 

Later being necessary for self aliened source ano drain 
structures. 

c) Another CVD process deposits an insulating layer (Si02) 
over the entire wafer, A maskinc step for defining areas 

of contact followed by chemical or plasma etching selectively 
exposes bare silicon or poly in the contact areas. 

d) Aluminium is deposited over the entire wafer by 
vacuume evaporation process. The next maskinc step patterns 
the A1 as desired for connections to electrodes. 



CHAPTER 3 


THRESHOLD VOLTAGE OF SK'ALL GEOMETRY MOSFETS 

3,1 INTRODUCTION 

Recent progress in lithography has led to the 
fabrication of MOSFETs with small channel dimension. As 
the device size shrinks, two departures from long channel 
behaviour occur. These departures can be studied indepen- 
dently and are referred to as the short channel and narrow 
width effects. Threshold voltage is found to show wide 
variation with device down scaling. To evaluate overall 
device behaviour and performance, it is necessary to model 
the threshold voltage to accomodate such higher order 
effects. 

In this chapter, first, effect of substrate bias on 
threshold voltage is considered. Two independent models 
have been presented for threshold voltage corresponding to 
short channel and narrow width effects. Finally all the 
three effects are combined to derive an expression for 
threshold voltage fox small geometry MOSFETs. The models 
are based upon simple analytical approximations. The 
emphasis is mainly on the simplicity, speed and' accuracy 
with a view to extending the MOS models for circuit simulator 
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To account the threshold voltage for graded substrate 
structure, a mathematical analysis has been presented to 
calculate the threshold shift* 

3*2 EFFECT OF SUBSTRATE BIAS 


The threshold voltage of a long channel MOSFET has 
been given in Chapter 2* 


^th ” ^FB + + - 




ox 


(3.1) 


Eq. (3,1) is threshold voltage expression without any 
terminal voltage applied and will have to be modified 
subsequently, 

When an inversion layer is formed in MOS structure 
there is a p-N junction between the surface inversion region 
(channel) and the electrical neutral bulk material. Effect of 
substrate bias is then to increase the depletion layer 
thickness in the bulk. 

Therefore depletion layer thickness on applying a 
substrate bias Vgg, 

= 

The threshold voltage expression (Eq, 3.1) Hx>difies to 


2 C-,^(20f^-VsB) 
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\h “ ^FB 


+ 20 ^ + 




ox 


A factor K may be employed to denote the effect of 
substrate bias 


'th 


= V, 


FB 


20 . 


2 €^c^^(20^) 


ox 


K. 


(3.2a) 


Vath Kg 


(1 + 

\ 



) 


3.3 SHORT CHANNEL EFFECT 


(3.2b) 


To achieve minimum feature length as the channel length 
is reduced, there is a considerable reduction in threshold 
voltage of MOSFET. The phenomenon, of short channel effect has 
been studied by a nunijer of researchers. Yau [7], Lee [8], 
and Merckel [9] have provided closed form expression for 
threshold voltage based upon two dimensional charge sharing 
theory. In this section a simple model based on the same two 
dimensional charge sharing approximation has been given. The 
model includes drain biasing effect. Schematic diagram of a 
short channel MOSFET is shown in Fig, 3.1, Departure from 
long channel behaviour can be obtained by applying charge 
conservation principle to the region bounded by the metal gate 
and the bulk of the semiconductor. 








20 




(3* 33 ) 


where, 

= total charge on the gate 

Q = total effective charge at Si-SiO« interface 
^ox 

0^ = total inversion layer charge "of mobile electrons 
Qg = total bulk charge due to ionized impurity atoms 
in the depletion region* 


Eq. (3.3a) may be expressed in terms of voltage, neglecting 
contribution due to at the onset of inversion. 





The threshold voltage is given by setting = 2 


^th = ^FB f 

ox 


(3.3b) 


with L = channel length 
W = channel width 

For short channel devices the full effect of on the 
threshold voltage is reduced, as shown in Fig, 3,1, This is 
because some of the field lines originating from the source 
and drain terminates at bulk. For small drain voltages 
(Vr,e < 10 mV), the depletion layer charge under the gate may 
be treated as contained in trapezoid AB’C’D. 
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The effective bulk charge is then 


Qo = KjQr 
^eff ^ ® 


where is the charge sharing factor 


Area (aB’C*D) 
“ AreaXABCD) 


Following the derivation given in Appendix A. 




L^2X. 


( 3 . 4 ) 


where 


^dm 


^ ^SB ^ 


c#l 








= VIIq, is obtained by solution of Poisson's equation, 
given in detail in Appendix A. 




+ (Wg+xJ^ ln( ) = 2 e . ” 

s j' r. ^ SI qN^ 


( 3 , 5 ) 
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For larger drain voltages, the depletion region near the 
drain expands further than that at the source. More charges 
under the gate now terminate at drain, as a result effective 
bulk charges under the gate decrease^ further. The charge 
sharing factor including the drain biasing effect modifies to 
(derivation in Appendix A): 


Kl = 


2L-[(rj+W^)^ - - 


^dm 




(3.6) 

Viith ’Wg given by equation (3,5), and given by; 


r,2-(r,+V)„)2 




The equation of threshold voltage 
MOSFET is modified as: 


qN^ 

(3.7) 

for short channel 


V+K = + 20^ + 


Kl«b 


where 


'th ~ ''FB ^ C~W 

o ^ 


Qg - qN^ W 
^fr' ^^s”^ 

qN ,X . L 

Vth = ^ =' Kl 

ox 


(3.9) 
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The factor being given by equation (3,6). 

In the above analysis, as the substrate bias is 
increased, the depth of the depletion region is also increased. 
Ref. to Fig, 3,1, point B* and C coincides together for large 
^SB* shape of the depletion region under the gate 

becomes triangular. 

The charge sharing factor in this case simply becomes 

ha If as 


_ Area (ABE) _ 1 
L ~ Area (aB CD) ~ 2 


where E is the point where B* and C* coincides. 
The depletion region width X^, now is 


2(L + 2rj) 

(3.10a) 

The expression for threshold voltage assumes the form 


Vth = VpB + 20p + 

ox 


(3.10b) 


Threshold voltage may be calculated using eq, (3,9) 
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•*8 2 

For <100> silicon substrate (Q = 1.4x10’’ coulomb/cm ) 
with aluminium gate (0^2 = - 0.92 volts) 

The flat band voltage 

Vpg = - 1.122 volts 

For long channel MOSFET (L=12 microns), with N^=2xlO^Vcc 

m * 

X . = 3.04 X lO”^ cms 

dm 

the threshold voltage 

Vth = 0.259 volts 

If Vgg = l.OV, 

= 0.674 volti^o 

The value of will decrease and even become negative 
for short channel devices. For circuit application, however 
a large value of is required. 

In the present case we have added a voltage term equal 

to 1.0 volt«. The origin of this term is related to ion- 

implantation with delta function approximation. The details 

* 

are discussed in Section 3,5. 

With this modification, theoretical plots of threshold 
voltage versus channel length have been given in Figs, 3.2 
to 3,4 , with drain voltage, substrate voltage and oxide 
thickness as parameters, 

A comparison between various models has' been given in 
Fig. 3.5. 



2. S0 



Fig. 3.2 
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CHANNEL LENGTH CMICRON£:> 


Fig. 3,3 


V 



2 



Fig, 3,4 


You Is m odel 
Mcrcket's model 
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CHANNEL LENGTH (MICRONS) 
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3.4 NARROW WIDTH EFFECT 

In contrast to the threshold voltage behaviour of 
short channel MOSFETs, narrow width -devices show an increase 
in threshold voltage. 

The cause of narrow width effect is chiefly attributed 
to the spreading of depletion layer charge under the thin 
oxide laterally [9]. With the reduction of the channel width 
of the MOSFET for VLSI applications, the effect assumes 
greater significance. An accurate expression for the threshold 
voltage should include process dependent considerations as 
well, which are described in this section. 

Fig, 3,6 shows the device cross section looking from 
the width direction. 

The gate overlaps the thick oxide on both the sides 
of thin gate oxide. The thick oxide on both sides does not 
change abruptly to the thin gate oxide, but is tapered and 
recessed. This results from fabrication process steps. 

The narrow width effect, which increases the threshold 
voltage as the width W is reduced, results from the charge 
stored under the thick oxide region and -the transition of the 
electrostatic potential from deep to the shallow depletion 
region. 
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As V< is reduced, the volume of the charge in the gate 
region is reduced, whereas in the tapered and thick oxide 
region the volume of the charge remains constant. This 
constant amount of charge becomes increasingly significant 
as V< is reduced, and contribute s to an increased threshold 
voltage. 

Rewriting the threshold voltage expression for long 
channel MOSFET, 

\h = VpB + ^ (3.11) 

here is the total depletion charge induced and is given by 

^%AP ‘‘‘ ^^th . 

where Qg = total charge under thin oxide 

= charge under tapered gate region 
Qth = charge under thick oxide region, 

Qg' is associated with a parallel combination of capacitances, 
the overall capacitance being 


C' = C’ + 2C-.- + 2 C,. 

T ox . TAP th 

Simple mathematical expressions are obtained for the various 
terms appeared above, the detailed derivations-are given in 
Appendix A, 
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qL X, M (b-a) 

Coulombs 

Qg = WL Coulombs 


Mm 


2 t3i(20f1-VsB) 


qN 


(3, 12a) 
(3.12b) 
(3.12c) 


C* 

ox 


ox 


WL 


ox 


(3.12d) 


where. 


^TAP 

- fd- in [ 

5 

2 d* b+f-, 
2 d'a+f*' 

. ( 3 . 12 e) 

d’ = 

1 +d^ 


( 3 . 13 a) 

H — 

^ox'*‘^res”^th 


( 3 . 13 b) 

U = 

(b-a) 


f 

2 dy' - 2 dg - 2 z' 


( 3 .i 3 c) 

g 

ftox+tres>a-tth-t 


( ’ 3 . 13 d) 


8 is the angle between capacitor plates given by 

^ = -tan-" C te)3 - tan-" [ 


(3.l3e) 



33 


ijue to the hioher threshold voltace, the depletion 
chsrce under thick oxide is sr.all cor.psred to that of 
and Qt-.-r, . Hence its con tricucion towards total depletion 

1 Mi" 

charge may be ignored without much loss of accuracy. 

The modified threshold voltace expression incorporating 
the narrow width effect, may be written as 


'^th ~ '■FE 


20, 


"ox'^^'^I.AP 



+ 20 £• + 



(3.14) 


where 






(3.1J) 


If a sucstrare bias is applied, the depletion region 
will spread into the substrate, but not under the thick ano 
tapered oxide region. Therefore lateral spreading width under 
these areas is approximately the width at the conduction at 
zero substrate biaso 

Threshold voltace calculations have been made using 
ec. (3,14) with following data. 
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t 

t 

t 


th 

res 

ox 

a 


b 

L 


= 1000 A® 

= 750 A® 

= 500 A® 

=2,0 microns 
= 2.5 microns 
= 25,0 microns 


Theoretical plots of versus channel width have been 

shown in Fig. (3.7), with doping density as parameter, 

3.5 SMALL GEa/.ETRY EFFECT 

The short channel effect and narrow width effect are 
in contrast to eachsother. To predict the threshold voltage 
for short channel and narrow width devices, both the effectr 
haj^to be combined, Yang [lO] presented a theory to calculat 
threshold voltage of such devices by treating the effects 
independently. 

For the present case, a general charge variation 
factor K is defined which is the combination of all the three 
effects described in this chapter. 


Thus 


'th 


= V, 


FB 


20 , 


2^se^A<20f) 

C 


K 


(3.16) 


ox 
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CHANNEL WIDTH CHICRONS5 


Fig. 3.7 




36 


Each factor contributing towards K, may be calculated 
as and when required, 

3.6 NON-UNIFORM DOPING 

In first order analysis of MOSFET, the doping concen- 
tration in the channel is assumed to be constant. However 
in practical devices, the doping is generally non-uniform. 
Even for substrates that are initially uniformly doped, 
because of thermal oxidation impurity redistribution, non- 
uniformity takes place. Above all in modern MOS technology 
ion-implantation is used extensively to improve device 
performance, which causes nonuniform doping near the surface 
region of the channel. 

In this section the effect of non-uniformly doped 
substrate on the threshold voltage is considered. The shift 
in the threshold voltage will affect the device and circuit 
performance. 

Substrate profile after Boron ion-implantation is 

shown in Fig. 3,8, The original implant alters, after 

thermal annealing. For a first order result of threshold 

voltage, the annealed profile is approximated by a step 

function, with a step depth X , equal to the sum of the 

s 



flunSTMArr nuClHli PUhPIIjK liKKP HUHUN JHI*l,*wrinij 

To't tMPI.AWTin,) 

STuP • 2 TTHr u II. I, 
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Fig, 3,8 Non-uniform Doping Profile 
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projected range and standard deviation of original profile 
(Fig. 3.9). 

In this case the inplant may be visualized as a delta 

function of negative charges (for Boron ion-implantation) 

localized at Si-Si02 interface. The charge is thus equivalent 

to a reduction in fixed oxide charges by an amount qD_, 

qE>I ^ 

v.'here D,. is the implanted dose, and this term = 1.0 volt 

O X 

has been taken in our calculations mentioned in Section 3.3. 

€30 

o 

N = Surface doping after annealing 
s 

N. = Substrate doping 


With the assumption that maximum depletion layer width 
under heavy inversion is less than X^, the surface region 
may be considered as a uniformly doped region with concen- 
tration N , the threshold voltage may be written as [11]. 

9 


^th 


'FB 


+ 20 , 


1 2 ^ qDj 


(3.17) 


ox 


ox 


In general the implanted dose of ions may lie partly 
in the oxide but the major part lies in the silicon. The 
portion of the dose in the depletion region alone is denoted 
as Dj 


0 




Fig. 3,9 Step Function Approximation 
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Relationship between surface band bending <f and gate 

s 

voltage '^GS may be obtained by solution of poisson' s equation 
[23], This becomes necessary specially when 


^dm ^ ^s » 

The relationship may be written as: 


^s exp(- ^ 


/Ys 


. {(^.m,.l)V2,^^2 3l/2 ( 3^,33 


where. 


a = 




ox 


,Lq 


Lq is debye length 

hi =1 


<^si^T 




y- = ^ volts. 

T q 


The parameters m^, m^ are the zero and first order 
moments of excess doping (N(x)--N^) of the non uniform surface 
defined as 


“o = fc / ^ (5^^ - 1) dx 


(3.19) 
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Xi 

- 1) X dx (3.19b) 

o A 


(where Xj is depth of non-uniform region such that N(Xj)=^^) 

A similar relationship for unimplanted or uniformly 
dpoed substrate may be obtained from (3,18) putting m^=mj^=0. 


V<3S„ = % * 


— [ (— 
Vt ^ \ 


exp(- 


-p ( ) 



(3.20) 


Under weak inversion condition, depletion layer 
charge per unit area for implanted subtrate is obtained as 

^ = |- [( -■>1-1)^''+ ^ ] (3.21a) 

OX ''t 

Similarly for uniformly doped case 

^ = |:- ( ^ -1)^/^ (3.21b) 

^ox 

Minority carrier charge at inversion may then be 
obtained by substracting depletion layer charge from total 
charge. 
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q ^Inv 

=ox 


= V, 


GS 




M( «yv2] 


(3.22) 


There is a considerable disagreement among various 

researchers about the threshold voltage definition and the 

surface inversion criteria for the non uniformly doped 

substrates [l2]-[l3]o For the present case however the threshold 

voltage is defined as the gate voltage required for which 

minority carrier density equals depletion layer charge per unit 

area devided by depletion region thickness. The threshold 

voltage shift is predicted as difference in gate voltage needed 

to maintain the same minority carrier density, N. in the 

inv 

implanted and uniform structure o 


Thus 


q N 


inv 


= V. - 


OX 


» f 


1/2 


m } 

# i 


<pg 

= V, - Cp - .a ^ 






- 1 ) 


1/2 


(3.23) 


Using above relationship, the threshold voltage shift 
may be written as 


A Vth = \ - \ 


'f’s 


s, Su'- ' V, 


1 _ .^l/2, “a Til 


= - % )+ H ^ -(^ -1) ] 


'fi 
“ 

T 

( 3 . 24 > 
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Fox a more realistic implanted profile, process simulator 
program SUPREM II [24] is used. It gives an accurate picture 
of the impurity concentration, obtained after a series of 
device fabrication steps. The threshold shift may be 

evaluated by calculating the parameters m^^ and m^. 

Considering a Gaussian profile implant of range R, 
straggle AFx, and dose D^, 


«(x) = 72^rW t 



- (R-t ) 

^ OX '' 

f2 AR 



Using equations (3,19a) and (3,19b) 



--[ 1+ erf ( 

Wk 


R-t 


ox 


f2AR 


)] 


(3.25a) 


m 


1 





. AR 


)(” ) 



exp[-(R-t^^)/(V2AR)^]^ I 
i+erfL(R-t^^)/V2^R ^ ) 


(3.25b) 


where erf(x) = ^ 



Surface potentials for implanted and uniformly doped 
substrate can be evaluated using following iteration scheme 
(derivation is in Appendix A). 
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I+l 


^ . 2 in . Xn [[ 14 ^ /(a 


.2[ ^!^/(a C„)] t( 4 i_„^.1)V2,^^J 


(3.26a) 


<p 1+'^ 

^SB 


N . 

+ 2 In ( -^ ) + In 
n . 




/(a C^,)]' 


q N. 

2 [ ... 


OX 


) 3 ( 




(3.26b) 

Threshold shift due to implantation may be evaluated 
using (3.24) with (3,25a), (3,25b), (3.26a) and (3.26b). 

As an illustrative example for a Gaussian implant with. 
Range R = 0.3 micron. 

Straggle AR = 0.08 microns 
Dose Dj = IxlO^Vsq.cm 

15 

is used for a substrate with = 5x10 /cc. 

Follovtfing results were obtained for a substrate bias of 
1.25 volts. 
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Threshold voltage shift 
A = 0<,415 volts 

A Wiide range of threshold voltage tailoring can thus 
he achieved by a choice of parameters R, An, and D^. 



CHAPTER 4 


operating voltage limitation due to scaling 


4ol introduction! 

A wide range of operating voltage is applicable for 
large geometry MOSFETs in terms of gate and drain voltages. 
However as the device size decreases, a number of physical 
mechanism tends to reduce the upper limit for operating 
voltage. These mechanisms originate from high electric field 
due to device dovwn scaling of oxide thickness and channel 
length. As the electric field is increased, the channel 
mobility becomes field dependent, when the field is increased 
further, carrier multiplication near the drain occurs leading 
to substrate current and parasitic bipolar transistor action. 
High field also causes injection of hot electrons into the 
oxide, leading to oxide charging and subsequent threshold 
voltage shift and transconductance degradation. 

In this chapter, some of the effects, limiting the 
operating voltage has been considered. These include: punch 
through effect, breakdown voltage, hot carrier injection, 
oxide breakdown and mobility degradation, A brief theory 
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has been given in each case followed by a simple expression 
fo- describing the behaviour of the respective mechanism for 
small geometry IwOSFET. 

4o2 Pin\ CH THROUGH EFFECT 

Among the aforesaid limiting factors for operating 
voltage of MOSFETs, punch through effect is significant as 
drain voltage onset of punch through goes down sharply with 
down scaling the channel lengths 

"Punch through* occurs when the drain space charge 
layer width reaches the space charge layer of the source. 
Upon occurence, shunt current flows in the bulk to augment 
the conduction current in the surface inversion region. This 
manifests itself as an increase in the output conductance of 
IqS vSo transfer characteristics [25]. 

An accurate expression for drain voltage onset of 
punch through may be obtained by the solution of Poisson* s 
equation in cylindrical coordinates. 

The Poisson* s equation in cylindrical coordinates 
can be written as 

6(r dv(r)/dr) q 



1 

r 




(4.1) 



Integrating with the boundary condition 


dV(r) 

dr 


l=(rj+ViQ) 


= 0 


or 


q q Ka 


(rj+W^)' 


2 € 


SI 


2 


integrating within limits (r^) to 


’'^bi 'hs = 2 6 ”. 

SI 


(4»2) 


) - 


^si 
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X In ( 






Vo = 


r^. ^ 5 ^ ^-Vfci-VSB 


(4.3) 

with Wg = L - V«g for punch through 

where the source space charge depletion layer thickness ¥*’ 
may be evaluated using eqr. » (3,5), 

As apparent from the expression itself to improve the 
punch through behaviour a high coping level is necessary but 
as this will affect breakdown voltage drastically, a trade off 
between punch through and breakdovn voltage is useful. To 



increase the punch through voltage single or double ion- 
ir.planted device structures can be made to increase the 
coping of the surface region, 

4.3 THE OXIDE 5REAKDa.»K' 

The dielectric strength of Si02 is of the order of 
70C7/p which lirrits the destructive gate breakdovi'n voltage 
to above 7 volts/lOC A° of cate oxide thickness, thus gate 
oxide thickness of the order of 200 A° are compatible Vvlth 
rr.oderate power supply voltages. However care has to be taken 
that the technological steps (in particular, the high 
temperature treatments) do not effect the Si-Si 02 (migration 
of the ions of the gate metal across the thin oxide layer). 

4.4 THE BREAKDOWig VOLTAGE 
# 

Apart from punch through, the breakdown phenomena 
also determines the highest applicable voltage, and limits 
the speed and power handling capacity of di| screte/integrated 
circuit device. Breakdown voltage is, therefore an 
important parameter in the design of the short channel MOS 
devices. 

The drain breakdovm voltage of a MOSFET is determined 
by the following two [ibj, phenomenon. 
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1* Avalanche breakdov.n of drain— substrate p-n junction 
caused by the hich electric field between the gate 
and the drain. 

2. characteristic with a negative resistance (switch 

back) originating from lateral bipolar action which is 
induced by weak avalanche hole current. 

The first one consists of simple p-n junction break- 
down theory and has been studied by several authors [loj. The 
modification due to lateral bipolar action is more appro- 
priate for small size (short channel) devices as will be 
apparent from further analysis. 

The hole current generated by weak avalanche multi- 
plication is injected into the substrate, so that some hole 
current can flow to the source. If drain voltage is love, 
most hole current flows out the substrate terminal. For large 
crain voltage, a substantial hole current can flov: to the 
source, and the product of the current and substrate resistanc 
becomes large enough to forward bias the source substrate 
junction causing electron injection into the substrate. This 
injection leads to a parasitic n-p-n ( sour ce- substrate-drain) 
bipolar transistor action (Fig. 4.1), The breakdown voltage 
of KIOSFET is then governed by the parasitic n-p-n transistor 



/ 


^c. No, 
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Fig. 4.1 parasitic Bipolar Transistor Action 
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hen floating base is assumed. From this point of view the 
cllowing bipolar equations can be used to explain the 
'.vitchinr phenomiena. 


ihe emitter (source) current 


I- = 


a D n -Y-r-/V-r. q D n -Voc/Vx 

= " -1) + A,s (e be' T, 

-(eff) "■ B 

( 4 . 4 ) 


ihc collector (drain) current 


q D n -V^c/V-r 
I, =. A„ -.-a-S (e BE' T _ 

i I 


.’here, 

A- 


""(eff) 




n 


V<- 

c 


= effective emitter area 
= emitter sidewall area 

= diffusion length of electrons in substrate 
= effective emitter length 
= the channel lencth L 


Effective current gain 
“eff - 


A 

L 'w_ 

n D 



(4.6) 


r ./L 


i nu s , 


a 


ef f 


TE 


ef f 


1-a 


eff 


L r . 

= — (-1 ) 


(4.7) 


The emitter collector breakoovm voltage of a npn transistor 
is given by following equation (with open base). 


BV 


BV 


CEO 


1 


fhe 


?E 


eff 


where BV^ is the pn junction breakdown voltage. 

Corresponding source drain breakdown voltage then 
Would be 


EVqs 






r . 


bo. (-J.) 


(4.8) 


JUNCTION CURVATURE EFFECT 


Drain breakdown voltage BVj is modified due to 
curvature of drain junction. Following [17], for cylindrical 
drain junction, the drain- substrate breakdown voltage is 
given by 
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Vcy = [| (Ti^+2Ti^/'^)ln(l+2n“®/‘^)-n^/'^]x BVj 

(4,9) 
r . 

where r] = ^ 
m 

‘V. depletion layer thickness at breakdowno 

As it can be seen from the eqn, (-^.S) the breakdown voltage 
is greatly reduced by decreasing the channel length L, This 
behaviour can be improved within a limit by using shallow 
souxce/drain diffusions and having larger values of L^. 

Theoretical plots of breakdown voltage versus channel 
have 

length £ been shown in Fig. 4,2, with the condition 
Lg = r j , and diffusion length ^ = 12o5p. 

4.5 HOT CAPIRIER EFFECT 

, Injection of hot carriers into the oxide causes 
performance variation during device operation. The physical 
mechanism is the emission of hot electrons from the silicon 
into silicon dioxide layer when applied voltage exceeds 
certain limits. 

Hot electrons can originate from the channel current. 
According to the ‘Lucky electron ivlodel’ [18], some of the 
electrons in the channel gain sufficient energy to surriount 
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Si-SiO^ energy barrier (barrier height c- 3.1 eV) without 
suffering an energy losing collision. These electrons are 
injected into the Si02» This phenomena is known as channel 
hot electron effect or CHEo 

The hot electrons can also be generated from substrate 
leakage current and is termed as substrate hot electron 
effect or SHE. 

Trapping of the injected hot electrons in the gate 
insulator layer can cause instability in the form of trans- 
conductance degracation and threshold shift. Threshold shift 
can also occur from hot hole injection but this usually 
occurs only when drain is heavily avalanching. 

In ordtrr to account for hot carrier injection and to 
understand guiding principles is optimizing device structure 
a deep insight into the mecnanism is necessary. A brief 
discussion on gate current and channel current will reveal the 
nature of hot carrier emission. 

The primary source of channel hot electron is channel 
current, as the voltage approaches drain breakdown voltage, 
^avalanche multiplication due to thermal leakage current near 
the drain also supply hot electrons in the channel. Electrons 
moving from source to drain gain energy from high field near 
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the drain and become ’heated electron gas'. A small fraction 
of such electrons surmount the Si-Si 02 energy barrier and 
are injected into the oxide. Under the condition V^-=Vr,c 
’lucky electrons' proves to be dominant than hot electrons 
caused oy avalanche multiplication. Thus the electron 
.rate current expected to be 



Is exp 


{- 



( 4 . 10 ) 


where , 

I = channel current 
s 

0^ = barrier height 

Tg = electron temp. 

The hole gate, current is generated from hot holes due 

to 'impact ionization' <, Most of the hole current, which are 
generated from avalanche multiplication near the drain 
under the bias flow towards the substrate as 

hole current. Only a small fraction of hole current is 
injected as hole gate current. 

Following the 'quasi thermal equilibrium approxi- 
mation' [I9j the gate current may simply be obtained by 
integrating hichardson's expression for the thermionic 


current 
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KT 


£ )V2 

27trr*^ 


exp 



(4.11) 


v;here , 

m* = effective mass of electron 
Hg = surface electron density- 

The electron temperature T^, can be obtained by using 
following expression [19J. 

OO 

"e = ( ^ (4.12) 

o ^ 5 T V 

e s 

where , 

Xg = electron relaxation time 

Ey = component of the electric field 

V„ = saturation velocity of electron, 
s 

An important assumption is that velocity V is parallel to 
the surface. 

For the small length devices the SHE voltage limits 
are much larger than CHE. As such the modelling of SHE is 
not much important, Hov/ever substrate current is of the form 

^sub ® (4.13) 


where b is a constant. 



Device Decradation: 


me hot carxiex injection into the gate oxice is known 
xox iir, posing, among the severest limitation on the miniatuxi— 
xation of MCSFtT in VLSI, This is because the part of the 
carriers trapped in the oxide cause long term device degra- 
dation. The initial damage appears a transconductance degra- 
dation, although very severe channel hot electron stress 
generally results in a combination of transconductance degra- 
dation end real threshold voltage shift. Long term operation 
of the device is seriously affected by oxide charging, because 
the charging continues to increase with time during device 
operation. As a result of this cumulative degradation, oxide 
charging limits the maximum voltage levels that can be applied 
fox a given specific device life time, 

iV.inimi ration of hot carrier effect: 

Fine line patterninc techniques such as electron 
beam lithography and reactive ion etching utilize radiative 
process steps. These process steps known to introduce charged 
as well as neutral electron traps in the oxide, which enhances 
device degradation due to hot carrier injection. Experiments 
shows that [32] threshold shifts due to channel hoc electron 
injection in devices fabricated by electron beam lithography 
is much more enhanced than device prepared by optical 
lithography. 
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A number of A'iOSFET structures were suggested as 'Hot 
carrier resistant' structures. These include As-P(n‘'' - n“) 
double diffused drain, drain offset gate and burried 
cnannel devices. The basic guiding principles for minimizing 
hot carrier generation are: 

1. Use of a graded drain junction for reducing electric 
field. 

2c Use of an offset gate for separating the gate electrode 
from the localized peak of the electric field, 

4.6 CARRIER MOBILITY 

One of the most important quantity which determines the 
electrical behaviour of MOSFET is carrier mobility. The 
carrier mobility varies significantly in small geometry MOSFETs 
due to higher internal electric fields existing in such 
devices. An exact model of carrier mobility is worthwhile 
since it determines the current levels in MOSFET, 

A brief theory for carrier mobility has been given 
including various influencing parameters on mobility. Later 
in this section a simple model describing exact mobility 
behaviour in MOSFET has been given. 

Mobility depends on factors such as the electric 
field strength, doping concentrations and the temperature. 



61 


However the variations dep€nd;> upon the band structure of 
the serrdconductor and the scattering mechanism which exists 
within a given device. The carrier mobility can be expressed 
by following simple relation 

q <T, > 

P = (4.14) 

where, 

'T = average carrier life time 
m* = effective mass of the carrier. 

The effective mass m* being a dependent on the band 
structure of the device. 


The impurity doping concentration and temperature 
influence on carrier mobility is ofcrvious as due to their 
impact on scattering processes. 

Few important scattering mechanisms are as follows: 


1, Coulomb Scattering: 


Occurs due to ionized impurities and most important 
of all processes at low velocities. Thus dominant at low 
temperature and high doping density mobility in this case [20] 


(m*)-l/2.T3/2 


N. 


(4,15) 


N 


I 


Ionized impurity density 
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2* Scatterinc due to Lattice Vibrationss 


This process is dominant at high temperature. The 
process induces optical and acoustic phenon scattering. For 
non-polar semiconductor like silicon the acoustic phenon 
scsxtering is more important y and the resulting mobility [21]: 


= 




(4.16) 


3, Ionized Impurity Scattering: 


Results due to electron-electron and hole~hole 
interaction. 


4o Neutral Impurity Cpllision: 

Results due to presence of neutral impurity atom 
having size different than host lattice atom causing a 
disturbance in periodic potential. 

5, D-fect Related Imperfections in Lattice. 

The other factors includes the radiation induced 
charges and interfacial states in the oxide degrading the 
surface mobility characteristics® 

All the abpve mentioned mechanisms contribute* 
together to produce overall mobility behaviour, however the 
exact value is determined by the operating temperature and 
impurity level* In ordinary conditions the net mobility is 
given by 
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+ ~ ) 




-1 


( 4 . 17 ) 


The electric field dependence comes about because of 
the chances in the momentum distribution among the carriers. 

This effectively transfers the carriers to different positions 
along the energy momentum curves thereby chancing their 
effective mass. Solving carrier velocity distribution function, 
srarting from the boltzman transport equation, foliovving 
relationship for the mobility results. 

q 

^ = 3 — ( 4 . 18 ) 

n:*<\r> 

.’..obility behaviour in MOSFET: 

As the current in MOSFET flows only in the inversion 
layer, mobility and drift velocity are expected to be influenced 
by the thickness of the inversion layer. From both field 
effect and hall effect it has been found that the free carrier 
mobility depends on the normal field i.e, the field perpen- 
dicular to the direction of current flow however it has been 

found experimentally [ll] that it is not a function of surface 

17 —3 

processing or doping density in the range I<;^< lo cm , At 
£ given temperature mobility decreases with increasing effective 
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transverse field, this observation is well supported by 
classical theory for surface scattering in a linear 
potential well. 

Another observed fact is at high lateral fields 

the rnobility also decreases due to velocity saturation, when 

the longitudinal field increases, eventually velocity 

saturation occurs similar to bulk silicon. For a given 

normal field E , the velocity is proportional to with 
A y 

a proportional constant equal to mobility. However as Ey 
increases, the velocity tends to saturate, 

A mobility model corresponding to the experimentally 
calculated values thus proposed as 

^ = o (4.19) 

( I+Sq 

v.'r.ere , 

= surface mobility under zero field 

0„ = Coeff, for normal field dependence 

u 

= Coeff, for lateral field dependence 
VjQ = threshold voltage with zero substrate voltage. 


For short channel MOSFETs the factors modify as follows: 
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a) Carrier Mobility under Zero Field: 

Baser on a nuiTibei of experinian tally verified 
results [22], the mobility under zero fields is modified 
s s 

^^o = "2 

where p- is bulk carrier mobility* 

b) The Coefficient G-:* 

Tnis parameter is a function of oxide thickness 
and is given as 


a . 

1 



where, 

a^ is a constant such that 

A Q 

a. = 10 cm /coul. fox p channel device 
1 ' 

7 0 

10 <a^<10 cm /coul. for n channel device* 


The Coefficient 6^: 

This parameter is a function of channel length and 
is given by a relation 

Sd = (L 

v.-'nere , 

E = Longitudinal critical field 



T«= 27 DEG C 


T* S® DEG C 
T» lee DEG C 


TOX=750A 
L=3.75E-4 
VDS=3 ,0V 
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For p channel device, ^ 10 V/^ 

For n channel devices < 3v/p 

Usin^ the relationships 




p = 


V, 


o 


( 4 , 20 ) 


Carrier mobility has a different temperature dependence at 

—2 4 

different ranges. Hov/ever for practical purposes a T 


variation can be assumed for n type silicon devices. 


Theoretical plots of mobility variation with gate 
have 

voltage [_ been given for three different temperature in 
Fig. 4.3, with following data: 

Pg = 1500 cm^/ VO Its- sec. 

a^ = 5x10 cm /coulo 

h = 2.5 volts/micron 
o 


and 


L = 3o75 microns 
V^S = volts 

= 750 A°. 

OX 



CHAPTER 5 


CO.>L>UCTia\' IK MOSFET AT BELOw Ai\'D ABOVE THRESHOLD 


IKTRODUCTIOK 

State of the art l\OS LSI, like random access memories 
(dynamic and static), microprocessors, read only memories and 
shift registers use I/.OSFET with channel length i to 6 microns* 
It is therefore, important to predict the conduction behaviour 
of MOSFET correctly for these range of dimensions* The binary 
logic states for the MOSFETs are defined as follows: 


LOGIC STATE 

DEVICE OPERATIQK 

NODE 

voltage 

Zero 

OFF 

^GS 

< Vth 

One 

OK 

^GS 

> ''th 


However there is a finite drain current even at < V.. . 

Via th 

Tiiis suDthreshold current or weak inversion current is 
important. In high performance dynamic RAMS, the magnitude 
of subthreshold current results in poor refresh times, 
pattern sensitivity and poor long cycle operation. Magnitude 
of subthreshold current becomes even more significant as 
channel length reduceso In the following sections a simple 
device model for subthreshold conductionis given for short 
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channel MOSFETs based on analytical approach and is more 
suited for computer aided LSI design^ 

The original long channel current voltage relationship 
has been modified for small geometry MOSFETs. All the small 
geometry effects discussed in the previous sections have 
been accomodated in the present model, 

5J2 SUBTHRESHOLD CQN’DUCTION 

For the gate voltage below threshold voltage, the 
semiconductor surface is in weak inversion. Corresponding 
drain current i.e, subthreshold current behaviour is, 
obviously totally different from that of above threshold. 
Subthreshold current for long channel MOSFETs has been well 
characterized by a number of researchers [26]-[28], but 
published work is limited for short channel MOSFETs. Literatu 
available for short channel case, mostly uses numerical 
methods and fitting parameter techniques [29]-[30]o These 
parameters often have to be re-evaluated for 

different devices or different ranges of values and therefore 
not suited well for circuit analysis or statistical modelling. 
A complete analytical model also provide better insight into 
the operation of device. In the present section one such 
analvtical model has been describedo 
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The drain current consists of two components 
nemely drift and diffusion current, 

^DS “ ^drift ^diff (--i) 

Idrift dominant only if the gate is biased above 
threshold and if 4 - 0, while is present with a 

carrier concentration gradient along the channel. The drain 
current in weak inversion is dominated by diffusion and is 
derived in the same way as the collectox current in bipolar 
transistor with a homogeneous base doping. 

Considering the MOSFET as an n-p-n (source-substrate- 
drain) bipolar transistor, 

IqS = a ^ (5o2) 

where , 

A = cross section of current flow 

D = diffusion constant of electrons in the channel 
n 

n = density of electrons in the channel. 

Since current is diffusion current and thus distribution of 
free carriers along the channel is given by the solution of 
following continuty equation: 
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An 


= 0 


(5.3) 


where. 


n^-n 


= diffusion length of electrons in the channel 

n being the electron density at source depletion 
s 

layer edge, 

;o> (b.3) may be solved using following boundary conditions. 


Anl = 0 

‘y<= 


(5.4) 


nli^y^ = Ter 


(5.5) 


where y and yr^ are the source and drain depletion layer 
S •lJ 


thicknes se s. 


From the solution of eq, (5,3), the free carrier 
density variation A n along the channel at the surface is 


given by 


An = ng(l-e 


y-y 

sin h( —T — ^ ) 

-qVps/kT ^ 

sin h( 


(5.6) 


electron density at source depletion layer edge, ng is 
given by 
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q %/^T 

® 

where , 

= electron density at gate depletion layer edge 

cp = surface potential 
f / kT 

Hp = n. e 
b 1 

The Surface Potential 

' 5 

Using the geometrical approximation for short 
channel devices as in Chapter 3, with refer to Fig. 3.1, 
the voltage applied at gate may be related to surface 
potential as follows: 

VqS = + V (2 631 q 

The geometrical charge sharing factor F has been obtained 
by same trapezoidal approximation, however, to obtain a 

simple close form expression fox V’ , the junction depth 

\!> 

Xj is assumed to be large, and F assumes the form 

(L - - Wj.) 

F = [ i + ^ 

where Wg and are source and drain depletion layer 
v^idth, and 
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qK 


''A 


(t.9a) 


Yn 



qKA 


(5.9b) 


i;ovv since <Pg, for subthreshold region varies over a small 

range ( 0.3 - 0.6), an average value of may be used 

in square root terms in eq. (5.9), due to its weak 

dependence on present case average value of 

is taken as = 0.45V. 

®av 

Using these assumptions, the surface potential 
becomes, (solving eqn» (5.7)), 

Cpg = -Vg- ) 


(5.10) 


with 


V‘ 


^si2V! 

““ c2 


OX 

with F given by eq, (5,8). 

Cross Section of Current Flow: 

Cross section of current flow A may 
A ■= W X 


be written 


as 
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where = channel thickness. 


depression for is derived in the appendix B, and is 


«si 




(^■.11) 


Subthreshold current is then determined at the edge of 
drain depletion region 


^DS “ 


y=L-YB 


Jsing equation (5o6) 


X,.n 

^DS = 


1 -qV-_/l^T 

s I (1 e ) 

* L-y -y ^ i 

tanh{ r ) 

Hi 

(5<.12) 


I!ov; for short channel devices 




since normally > 10 


drain subthreshold current is then 


( L-Yj-YoT 


•q 0f/KT -qVos/KT 


[1-e • ]x expt-jjP ] 


(5.13) 
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■^DS 


vs. 


'DS 


characteristics - 


i) Fox small drain voltace ( ^ few ~ ) the value of 

q 

is taken as constant since dependence of V-.^ on ^ comes 
only through the charge sharing factor F. The varies 
exponentially with 

ii) For voltages > 3, ~ , the exponential term in 

Vjs» rapidly vanishes, and variation of y^ and dominates 
Thus in this range, leaving the exponential term in 


I 


DS “ 


q’' °n ^1 
( L-Ys-Yd) 


-q0f/KT 


q 

® [ "kt” ^ 


(5.14) 


Eq, (^,14) shows cross product of drain and substrate 
voltage, hence as the substrate bias is increased, the 
dependence of current upon drain voltage will also increase. 
However, the dependence of the current on drain voltage 
will be a function of oxide thickness and- substrate doping, 

V«ith increase in drain voltage VuS- drain depletion 
layer width y^^ increases and as predicted by eq, (5.14) 
drain current also increases for short channel devices 
the increase in drain current is much more prominent, and if 
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if the channel is sufficiently short, that for a certain 
ora in voltage 


i-^Ys + y^, 

condition holds v^ithout occurence of avalanche breakdown 
in drain, the drain current is wholly governed by denornenatc 
term in eqn. (5ol4) that anomalously large current is 
predicted at punch through. 


This an maly results due to dominance of drift 
component of current at drain depletion layer edge, which 
has been taken to be negligible in earlier assumption 
through the boundary conditions. Assuming the transition 
from diffusion to drift occurs./ over one debye length, 
the drain current at punch through 


I 


D 


qD W K e 
^ n 1 


-q0f /KT 




exp [q^g/KT], 


/KT 

q ^ e exp[q ] 


I i 

'fs) 1 

2tsi<''bi-V''Ds) , 1 



<c 

cr 

i 



( 5 . 15 ) 
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At the punch through, punch through voltage may be 
determined by using expression given in Chapter 4, 

Theoretical plots of subthreshold current versus drain 
voltage have been shovm for fixed gate voltage V^r-=0.2 volts. 
The long channel threshold voltage being 0.5 volts, 

5.3 above threshold CaxDUCTIO; 


Current voltage relationship at saturation and linear 
region of operation has been discussed in Chapter 2, In fact 
the manner in which the MOSFET has developed, the observed 
short channel effects are not explicable from the normal 
gradual channel approximation. In addition to various 
consequences of twr© dimensional sharing of substrate 
depletion charge between the gate, source and drain regions, 
the charge sharing affects above threshold behaviour of 
LCSFET both in the linear and saturation conduction 
characteristics. 


To clarify the situation, normalized drain current 

I'^ is defined as 
u 



■DS W 


(5.16) 


where L and 'w are the effective channel length and width 
respectively. As long as the conventional long channel 
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theory holds, the vs. do nor depenc on L and Thus 

“ '^DS 1® required to be mocified to include 

vc.rious scaled dov-.n effect,, 


The overall current dependence on the drain and gate 
voltage has been taken to be of the followino nature. 


i) At linear or triode region of operation 


"DS 


ox 


V? 


fCv „-V* )v - 1 


(5.17a) 




ii) At saturation 


C W 
ox 


^us = i'* aTETSI) t'^GS - 


(5,17t) 


'^us ' ''gs - V|h > 0- 


Where p,* effective mobility discussed in Chapter 4, 


Vf, = effective threshold voltage oiven as 
th " ' 


VJh = Vr;, + 20f + 


2t^iqig(2 )6,)'.K 


TB 


+ AV, 


ox 


th 


(5.18) 


The factor K consists of one or more factors corresponding 
to various scaled down effects such as short channel effect, 
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narrow width effect anc. substrate biasing effect, as 
dfc. scribed in Chapter 3. 

Vth induces the shift in threshold voltage, due to 
non-unif orr.ly doped substreteo 

Finally the rerri. AL, known as channel reduction term 
is determined as follows. 

The Channel i^eduction Term: 

tihile operating in the saturation region, the output 

characteristics of J%OS devices shows a finite differential 

on-resistance, the potential at the end of inversion layer 

will be fixed at the value Vp. o As the drain voltage is 

^sat 

increased further, the reverse bias across the drain junction 
'.viil increase, as a result the depletion region separating 
end of channel and drain will increase. Effective channel 
length thus will become shorter, resulting in an upward 
riltinc of the current voltage characteristics beyond 
saturation. As it is obvious from the discussion, this 
phenomena is particularly important for short channel device. 

Out of the various models given for predicting the 
channel length modulation term [31], at present an analytical 
raodel based on velocity saturation effect has been chosen. 
This model is more suitable for short channel MOoFEi s. 
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This mooel is based upon the assumption that carrier- 
vclocitifes saturate in the channel* The velocity saturation is 
apparently, field depenoent and occurs at some critical field 

strencth L . 

0 

For simplicity it has been assumed that the channel 
of 1:03? LT in saturation jMODH is composed of only two recions. 
The normal channel, and the space charge region L, as shown 
in the Fig* to2. 

In the space charge region there is a two dimensional 
current flow near the drain and carriers are spread in the 
region POK. This bulk conduction effect may be attributed 
to two dimensional electric field distribution. As a 
consequence of this bulk conduction, it becomes necessary to 
include the contribution of the mobile charges to the drain 
substrate depletion charge. Including it in Poisson’s 
equa tion 


£ilxi = _i 

cly2 


^si ^ '^sat 


where , 

J(y) = current density 
V . = carrier saturation velocity 
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the differential equation may be written as 


d^V(y) 1 , .. 

^ = t — (q.^. + 


set 


dy 


-SI 


Vi/ •V .X . 

sat j 


(^•.19) 


and may be solveo using following boundary conditions 


gy 

dy 


L~ AL 


= E(y)l 


= E 


L- AL 


V(U =Vds + v^, 


V(L- AL) = V, 


D 


sat 


Thus V(L) - V(L- AL) = (V ^3 + ) 


sat 


Inteoration of eqn. (5ol9) twice with boundary conditions 
result in a quadratic equation in AL. 


1 f X, ^sat 

^si ^ ^ ^'"•^j-'^sat 


) (AL)^ + 2^(^U+(Vq3-kV^.-Vq ) = O 


sat 


(5.20) 


A L = 


L may be obtained by solving (5.20) 

I 

i- (c 3 N.+ 

^si 




D. 


1 , ^ sat 

€si A*^ 


) 


(5.21) 
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The- current I fronr. (v.l7b) 

"^sat 


(v,s-v|h )" 


(5.22) 


A si^;ul^aneous solutions of above two equations will 

et ‘■rrr.ine the saturation current level l_ and channel 

"^sat 

length reduction A Lo 

Ij^S “ VqS characteristics have been shown in Figs. 

(5o 3) - (5.4), based on computer calculations. 

The effect of channel length reduction on current voltage 
characteristics has been shov-Ti in Fig. 5.3. 

Fig. 5,4 shows the effect of channel width reduction 
on current level. 


For both the cases the ii/L ratio is fixed. 



i o 


'LJ: 


;c. OF HI.- 


u.»Oiir. ErrLwTi \j.\ ulrlCUIT PERFOruVlAI'JCt 


o • X Xi't i i iC'ijO «,/ i Z 0, H 

Continuin:. irriprovement in ?.1CS fabrication technology 
have made possible steady reduction in the internal aimen- 
sions of semiconductor devices. The trend towards progre- 
ssively smaller internal dimensions is likely to continue 
for some years until fundamental physical limitation of 
transistors are approached. As the device size reduces, 
it affects the various parameters governing the circuit 
performance of the device. In this chapter influence of few 
such parameters on circuit performance have been considered. 

The basic circuit chosen for the discussion is the 
simple MOS inverter. It is sufficient to serve the present 
purpose, since tt\e inverter is the basic circuit from which 
most MOS logic circuits are developed.. The MOS inverter 
exhibits all the essential features of MOS logic gates 
except for logic functions. Subsequent extention of MOS 
inverter concepts to KOh and cate is simple. 
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6.2 INFLUENCE OF THRESHOLD VOLTAGE 

The thresholi voltage behaviour for short channel 
jV.OSFET has been discussed in Chapter 3 in detail. Apart 
frorr. the wide variation of threshold voltage with device 
down scalin-:", it also varies considerably with operating 
point of thf circuit due to drain and substrate biasing 
effects. In this section, effect of threshold voltage 
variation due to device ninimi ration effect on inverter 
characteristic has been given. 

The load device, in an inverter is used as a current 
source. Its threshold voltage is tailored by ion-implantation 
to provide a specific saturation current level. For simpli- 
city this device is assumed to have negligible effect on 
overall characteristic due to threshold voltage shift. Impor- 
tant parameters such as noise margin, logic threshold and . 
logic swing are then governed by threshold voltage of driver 
device. A variation in threshold voltage of driver device 
will reflect its effect on the entire characteristic of the 
inverter. 

As an illustrative example, two inverters have been 
Considered with different channel geometries. The depletion 
load device has been taken same for both cases also the 
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acp^'Ct ratio (n/L) is same for both the inverters. The 
parameters for io-i'. device are as followsi 

Jnennel length L =20 microns 
Channel v.iath ki = 5 microns 
Uniform substrate doping K. = IxlO^^/cc 
Gate oxide thickness T ^ = 500 A° 

O X 

Thresholo voltage = -2.5 volts 

The parameters for driver device are listed in Table 6.1. 

Threshold voltages have been calculated by using 
models developed in Chapter 3„ Transfer characteristics 
of both inverters are drawn analytically and plotted in 
rig, 6,1, 

Calculated specifications for both the inverters 
are listed in Table 6.2, 

Apparently there is a considerable change in 
inverter parameters due to down scaling. 

6.3 INFiajB'JCE OF SUBTHRESHOLD CURRB'IT 

Due to subthreshold conduction mechanism, MOS 
transistor drain current does, not go to zero at 
but falls of exponentially. For long channel devices, 
this current is, of the order of a few tens of pico amps, 



Vour — »• 



Fig. 6.1 Inverter 


Transfer Chsrecteristics 


Table 6.1 


Lasic set of specif ication used for the Illustrative Example 


larametrx 

1 Large geometry 
devi ce 

II Small geometry 
device 

Channel length L 

8 microns 

1 micron 

Channel wicth «• 

64 microns 

8 microns 

Uniform substrate 

doping 

IxlO^Vcc 

1x10^^/ cc 

Gate oxide thickness T „ 

OX 

600 A® 

500 A° 

Junction depth r. 

0.5 micron 

0.5 micron 

Substrate bias 

4 

l.OV 

l.OV 


Table 6.2 

Calculated specifications for inverters I and II 


Specification 

Inverter I 

Inverter II 

Threshold voltage of 
Driver device 

Oc96V 

0.59V 

Logic threshold 

1.4019V 

1.03194V 

Uoise margin 

1.15 V 

: 0.6V 

Uoise margin NM2 

3,4V 

3.8V 
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hc-nc€ insi ;inif For short channel length devices, 

cubthies- old cuti^nt rices to much higher level, even to 
the ordcx o: microemps. The effect of subthreshold current 
is more serious for cynamic circuits, where the magnitude 
of subtnresi.old conduction current results in a poor refresh 


times anu poor long cycle opera tion*, For a reasonably safe 

operation, thf value of subthreshold current at =0 is 

-5 

dtsioneo to be at least 10 times smaller than the current 


around Vq^ == 

Among the various other higher order effects, the 
field dependent mobility phenomena plays a vital role in 
circuit performance. The reduction in current levels due 
to this effect increases the noise margin, while the 
capacitive charging and discharging time becomes higher. 



CHAPTER 7 


s V Uj oi el'll 

A cad niodC'l of a I'.OS transistor is a trade-off between 
sirr;plicity and accuracy. For earlier lone channel devices, 
the gradual channel approximation with one dimensional theory 
v.as satisfactory. However with the advancement in technology, 
devices have become smaller and a two dimensional analysis 
with other higher order effects have become necessary. In 
fact only a numerical solution method achieves sufficiently 
hi: h accuracy but then at the cost of simplicity. On the 
other hand, integrated circuit designers continuously ask for. 
simple but accurate device models in order to have aids 
amenable for circuit simulators. 

The latter requirement has been the motivation for 
present v*ork, where we have tried to focus our attention in 
studying and understanding various physical phenomena, 
originated as a result of scaling dovm the device size and 
also development of several proposed analytical and pseudo- 
analytical small dimension MOSFET models. To minimize the 
undesirable impacts of higher order effects a careful and 
proper selection of device parameters is required, A deep 
understanding of physical aspect helps to know each such 
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ptir3rnPl«©x sno dlso idsntifiss th© various fsctoxs which 
control tht /..CSFilT electrical behaviour. These parameters 
v.i'wh the accurate moocl ct.n be usc-a for integrated circuit 
ccsign. 

The models which wo have developed for the scaled 
Qov.TJ versions of the K'uSFETs, therefore, assume sufficient 
iraportance. They link various electrical characteristics 
of such devices, directly with physical phenomena and of 
course a better insight into the mechanism is an added 
advantage. For example, the theoretical plots in Chapter 3, 
give a first order estimate of an important parameter, the 
threshold voltage for various combinations of other 
parameters, valid upto very small device dimensions. Simi- 
larly the models in Chapter 4, provide a guideline for 
safe operation of the device. 

We can suimnerize the important points of the 
Work into following steps: 

1, For the short channel MOSFETs, the threshold voltage 

decreases as channel length decreases. For such devices, 
threshold voltage decreases as substrate ooping concen- 
tration, gate oxide thickness and substrate bias voltage 
but increases v;ith decreasing drain voltage and junction 


d ep th , 



2* Fox the narrow width MOSFETs, the threshold voltage 
increases as channel width decreases. This influence 
is higher at increased substrate doping concentrations. 
Hence the threshold voltage sensitivity to the substrate 
bias is increased for such devices. 

3o For small geometry MOSFETs punch through occurs more 

easily as the substrate doping concentration and channel 
length decreases. The required drain voltage values for 
the punch through are larger than the required substrate 
bia s. 

4o The source drain breakdown voltage of the short channel 

iViOSFET is greatly reduced by decreasing the channel length 
Lo To improve the bxeakdov^mi voltage the resistance of the 
substrate can be minimized, so that potential drop 

across K . due to substrate current remains smaller than 

SUD 

C.6V, The use of snaliower source and drain reduces the 
probability of corner breakdown phenomenon. 

5o To reduce oxide charging due to hot carrier effects, the 
density of vaster related traps in the oxide should be 
minimized, because such traps are knovm to capture 
electrons. Design of a circuit consistent with the hot 
electron limitation is for same circuits, nontrivial. 

The 1 pm technology may be expected to run with a drain 



supply voltage of about 2.5V. Foi conventional 
circuits this voltage sets the upper bound on both 
and Vpg. Thus circuits are unconditionally safe as hot 
electrons are concerned. Also extremely shallow source/ 
drain junction is found to increase hot carrier junction 
phenomenon. 

6. For short channel devices the subthreshold current depends 

on drain voltage fox VQ2,>3KT/q, unlike the lone channel 

devices. This may be attributed to early punch through 

in such devices. An appropriate increase in doping level 

thus, helps checking large subthreshold current to an 

extent. The subthresr.old characteristics are important 

for utilizing devices in logic circuitry. However simple 

is 

linear scaling of a device /_in sufficient to hole subthre- 
shold current to an acceptable level, K'any factors such 
as band gap, temperature, interface charce-and diffusion 
profiles are not scalable parameters but they have a 
profound influence on device performance in subthreshold 
region. 


Finally, we conclude with a mention of the scope of 


future Work on this tonic 


The aspects which we have covered 



here, appropriate numerical models can be developed for them, 
based upon the study already made^ Such models will give 
more accurate results at a cost of increased computation time, 
k closed form relationship for hot electron effect may also 
be useful. Another useful and interesting area to work on 
will be development of computer programs for the transient/ 
small signal analysis and automated parameter selection of 
scaled down MOSFETSc 


I 



appendix a 


Short Channel MOSFET Calculations: 


The relevant diagram is shov.Ti in Fig. 1 , The 
charge sharing factor F is given by 

F - Area( aB* C* D) 

“ AreaCAbCDl 


Lgff + L‘ 


2L 


ef f 


From simple geometry 
■.«e get for small drain 


4 = 

L’ = ” 2L 


voltage (L^ = L2) 


- « r, - 

am-* 3 s 


( 1 ) 


hence , 


r 




am 


L-2X. 




For a larger drain voltage L2 > 
the charge sharing factor 


( 2 ) 
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APPENDIX A 

Short Channel MOSFET Calculations: 

The relevant diagram is shov.Ti in Fig, j 
charge sharing factor F is given by 


F _ Area(AB«C«D) 
“ Area (AS CD'} 


Lgff + L' 


2L 


'ef f 


From simple geometry 

Vie get for small drain voltage (L^ = 1^2^ 

L, = [(r..+Vi - r . - X, 

1 3 s' am"* 1 5 


L’ = Lg.f - 2L 


hence , 


r = 


[(r..+VV - X,^ 3^/2 

^ j sj am ■* 


L-2X, 


- r.-X 
L-_s. 


For a larger drain voltage > Lj^ 

the charge sharing factor 




. The 


( 1 ) 


F 


1 


(3) 




IOC 


i-.r-ain froir. geometrical considerations 

+ X|,_ = (‘»®) 

(Xjj+r^+Lj)^ + = (rj+'Vjj)^ (4b) 

Solving equation (4) 

h = [(=^3+Ws>^ - 

4 = [(rj+Wn)2 . 


The charce sharino factor F becomes 


F = 1 . ^ 4.n3 r^-s-^ 


D 


2( L-Xg-X^) 


( 6 ) 


If the backcate bias is hi: h enough the surface of 
the depletion region under the cate electrode become 
triangular* The depth of the depletion region at the onset 
of this transition (ioeo L'=0) can be determined by 


Li + L2 - 


( 7 ) 


Substituting equation (f) in (7) and solving for 
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2(L+2r .) 

* J 

( 8 ) 


Derivation of ’»*_ and Vv.: 

b xJ 


In cylindrical coordinates Poisson’s equation is 
v.-ritten as 


1 

r 






dr 


SI 


i.e ^ 


dr 




2 ^ 


+ c, 


Si 


v.'ith boundary condition at the source end 


Mil 

dr 

CO 

II 

H 

= E(r^+'Ag) = 0 

dV(rl 


(r +V. 

^ -IS 

dr 

■ 2^31 

” 2 £ . r 

^S1 


(9) 


Intecratinq (9) between limits r. to r 

J , j s 


^ si tVbi^VsB) _ 

qN^ - 2 + 


r .+V«' 

ln(-X_^) 


r . 
J 


1/2 


( 10 ) 
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Similarly for drain side 


€ . (V, .+Ven+V~e) 
^ SI bi SB DS*^ 




O r^ + n. 

+ (rj+V,-)2 In 


( 11 ) 

Source and drain depletion layer width V« and Vv~ 

o LJ 

may be evaluated from equations (IC) and (11) by any 
suitable iteration method. 


Narrow Vdidth Calculations: 

The schematic di'agram is shown in Fig. 2. The tapered 
oxide region extends over a distance of (b-a) units in the 
horizontal direction v;hile the Si-Si 02 place has been taken 


at a height of Y. 


res 


The charge contained in the triangle is 

expressed by following integral. 


L b 


Tap 


= q f S I ciy dz 

o a y^ ^ 


''A ""dm ^ 


couls 


( 12 ) 


Charge under thin oxide is as usual 




( 13 ) 
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CjAP calculated by considering capacitor 

model shovjn in Fig^ 3, The capacitor is 3 series 

combination of infinitesimal small capacitances as shov-n 
in Fig. These capacitances are further devideci into 
infinitesimal small parallel capacitances; 

One such combination of parallel capacitance may 
be written as 


dc = / 

x=R-t 


x.d£ 


(14) 


The total capacitance is thus 


Evaluating R, t, and b geometrically 


(15) 


r - U"d' r 2ti ’b+f i 

‘^AP “ ^ox “7 C 2d'a+fJ 


C16) 


here d' 


(17a) 


^Qx'^^re s'~^ th 
(b-3) 


(17b) 


= 2dy' + 2dg - 2z* 


(17c) 


(t +t ).a-t.. ob 
lox_resl - th 


(b-a 


(17d) 




where coordinates y* and z’ are given by 


y 


z* 


ireslith 

Hh-^ox 


^ox-^-Sh-° 


5 , the angle between capacitor plates is 

5 = _ tan"^ (^fp) - tan"^ ( %/ ) Cl7e) 

Derivation of Surface Potential; 

From eq. (3,23) 


ql\ 


inv 


'fc- 


= V- - ( 

C ^ ''Gt '^S- Vt ^ Vt 

ox I 1 1 , T 


® ' ' ^1 _ 1 \ 1/ 2 o 


- m^-l) 


IB 

+ ^ 


(IS) 


and 


2^nv ,, vA a 

C V ” ” S,, V_ 

o X 1 






v_ 

i 


- 1) 


( 19 ) 


Eliminatina from (IE) and (3, IS) 


qN._„ n. o VeQ 

^ = a [(.^ F exp (- ) exp(-^ 

T ox A T i 



( 20 ) 


~ nij- 1)^/2 + mj\^2 3^/^ 


< S / 

" ^ -^l-l)^'^+rr.^/'f2 1 


Similarly eliminating V- from (19) and (3,20) 

T 


u 


^ (- p>=p( 5 1/' 

'ox* '^T '- ^A ^T ^T ^T -* 
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Taking last term in (20) to left, squaring and taking log 
an iteration scheme fox finding results for (I+l)th 
iteration. 
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Simi larly, 
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'T ox 

A few iterations using (22) and (23) gives and^- result: 

V7ith high accuracy. 





appendix b 

Calculation of Channel Thickness for Subthreshold Conduction; 

To find the channel thickness X^, it is necessary to 
consider the total charge available for conductibn in the 
surface depletion region. By integration in the X direction 
over the surface deplexion region the total charge is 
obtained as 



The free charge available for conduction is thus 
confined to a distance from the surface of 


X 



v.’hich is the channel thickness 
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